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1 Introduction 



In the revolutionary braneworld viewpoint, our universe is a 3-brane embedded in an extra 
dimensional bulk. Standard matter and all interactions are confined on the brane; only 
graviton and possibly non-standard matter are free to probe the full bulk [1,2]. Based on 
the braneworld viewpoint, our universe may contain many more dimensions than those we 
experience with our senses. The most compelling reasons to believe in extra dimensions are 
that they permit new connections between physical properties of the observed universe and 
suggest the possibihty for explaining some of its more mysterious features. Extra dimensions 
can have novel implications for the world we see, and they can explain phenomena that seem 
to be mysterious when viewed from the perspective of a three-dimensional observer. Even 
if one is doubtful about string theory due to, for instance, its huge number of landscapes, 
recent researches have provided perhaps the most compelling argument in the favor of extra 
dimensions: a universe with extra dimensions might contain clues to physics puzzles that have 
no convincing solutions without them. This reason alone makes extra dimensional theories 
worthy of investigation. In this streamline, the braneworld models that are inspired by ideas 
from string theory provide a rich and interesting phenomenology, where higher-dimensional 
gravity effects in the early and late universe can be explored, and predictions can be made in 
comparison with high-precision cosmological data. Even for the simplest models of RS and 
DGP, braneworld cosmology brings new implications on the infiation and structure formation 
[1,3]. Also it brings new ideas for dark energy and opens up exciting prospects for subjecting 
M-theory ideas to the increasingly stringent tests provided by high-precision astronomical 
observations [3] . At the same time, braneworld models provide a rich playground for probing 
the geometry and dynamics of the gravitational field and its interaction with matter [3]. In 
these respects, the braneworld model of Dvali, Gabadadze and Porrati (DGP) is a scenario 
that gravity is altered at immense distances by the excruciatingly slow leakage of gravity off 
our 3-brane universe. In this braneworld scenario, the bulk is considered as empty except 
for a cosmological constant and the matter fields on the brane arc considered as responsible 
for the evolution on the brane [4,5]. The self-accelerating DGP branch explains late-time 
speed-up by itself, without recourse to dark energy or other mysterious components [5,6]. 
Even the normal DGP branch has the potential to realize an effective phantom phase via 
dynamical screening of the brane cosmological constant [7]. 

Here we are going to study cosmological dynamics in a DGP setup with a bulk canon- 
ical/phantom scalar field. Many authors have studied the cosmological consequences of a 
bulk scalar field (see for instance [8,9,10]). One of the first motivations to introduce a bulk 
scalar field was to stabilize the distance between the two branes in the Randall-Sundrum 
two-brane model [11]. A second motivation for studying scalar fields in the bulk is that such 
a setup could provide some clue to solve the cosmological constant problem [12]. Models 
with infiation driven by bulk scalar field have been studied and it is shown that infiation 
is possible without infiaton on the brane [13]. Generally, the scalar field living in the bulk 
affects the cosmological dynamics on the brane considerably. The evolution of this field has 
some interesting cosmological implications; it can give rise self-acceleration and phantom-like 
phase even in the normal DGP branch of the model in some appropriate situations. Solving 
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the field equations for a braneworld cosmology with bulk scalar field is not generally an 
easy task. Nevertheless, during the past decade attempts have been performed to handle 
this problem. As an attempt, one can express the five-dimensional Einstein equations in 
terms of 4-dimensional tensors on the brane. Then, by using the Darmois-Israel matching 
conditions, one can determine these tensors [14]. However, this approach cannot determine 
all the tensors on the brane. As has been pointed out in Ref. [15], one can proceed further 
by making assumptions about the bulk solution. But, it is not obvious that such assump- 
tions are justified [9]. To determine which solutions of the 4-dimensional Einstein equations 
are allowed, one should solve the full 5-dimensional field equations. Full bulk solutions for 
static braneworlds have been found in some special cases (see [9] and references therein). If 
the scalar field potential takes a supergravity-like form, the field equations can be reduced 
to first order equations, which can be solved with relative case [16] (sec also [9]). In this 
paper, we generalize the work of Davis [9] to the DGP setup. We consider an extension of 
the DGP scenario that the bulk is non-empty and contains a canonical or phantom scalar 
field. Since the self-accelerating DGP branch has ghost instabilities, we restrict our study to 
the normal DGP branch of the model. The bulk equations of motion are derived and some 
special classes of solutions are presented. We determine the evolution of the brane when 
the potential of the scalar field takes a supergravity-likc form. Some clarifying examples 
along with numerical analysis of the model parameter space are presented in each case. The 
importance of this work hes in the fact that bulk scalar field in DGP setup has not been 
studied in supergravity-style analysis yet. Also our detailed study in this framework reveals 
some yet unexplored aspects of cosmological dynamics of the bulk scalar field in DGP setup. 



2 A canonical bulk scalar field in the DGP setup 
2.1 The bulk field equations 

The five-dimensional action for a DGP-inspired braneworld model with a bulk canonical 
scalar field can be written as follows 

[ d'xV^{^('^R-lm'-V{cf>)}- [ (i^xv^(-^i?+-^[X]+A(0)) , (1) 

J bulk -^l^t ^ J brane ^ 2.K% ZK.^ ' 

where qab is the bulk metric and Uab is the induced metric on the brane. They are related 
by Uab = Qab — fiAUB, where ua is the unit vector normal to the 3-branc and A, B are 
the five dimensional indices. The Gibbons-Hawking boundary term is included via jump of 
the trace of the extrinsic curvature [K] in the brane action. Also, k-I = where M5 is 
the fundamental five-dimensional Planck mass. The brane Lagrangian Cb{(f)) includes all the 
Standard Model fields which are confined to the brane, and depends on the bulk scalar field. 
Varying the action with respect to the bulk scalar field and also the bulk metric gives the 
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bulk equations of motion 



.2, dV V^dCt{(l)) 



where 



rp(brane) r) ^-^brane , r, r ( a\ 

J-AB — ~^ ^^AB + f^ABi-brane- {^) 

Ttr"^ is the energy-momentum tensor locahzed on the brane. is the Dirac delta 

function with support on the brane which we assume to be located at y = where y is the 
coordinate of the extra dimension. The action (1) implies the following jump conditions 

N Va(?\ = , (5) 

KAB-KhAB\^->^lTTB''''\ (6) 

The Gauss-Codacci equations relate projections of the bulk Einstein tensor to the extrinsic 
curvature via 

N''GABh^ = DAK^~DcK, (7) 

where Da is the covariant derivative with respect to the bulk metric. Combining this with 
(3) and substituting it into the jump conditions (5) and (6) we find the effective energy- 
momentum conservation equation on the brane 

To formulate cosmological dynamics on the brane, we assume the following line element 

ds^ = gABdx^dx^ = -n^{y, t)df + a^{y, t)-fijdx'dx^ + b'^{y, t)dy^ (9) 

where 7^ is a maximally symmetric 3-dimensional metric defined as = 5ij-\-k-^^^ where 
k = —1, 0, +1 parameterizes the spatial curvature and = a:^a:^ 

Since we consider here homogeneous and isotropic geometries inside the brane, 
can be expressed quite generally in the following form 

rpA{brane) ^ ^^■^g(^_p^ , , Pb , Pb , 0) . (10) 

The extrinsic curvature tensor in the background metric (9) is given by 
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So the jump conditions (5) and (6) are given as follows (we note that in the forthcoming 
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equations, ^ = Tc where Tc is the DGP crossover scale) 



\a'] nl^ , _^/ao , u^l 



4 = t(3P' + 2/'.) + -?.(-4-2^ + ?!^-4| (13) 
noOo 3 kIuq I Oq aono ao Oq J 



[01 _ 5£6(0) 



(14) 



6o (^0 

where a prime marks differentiation with respect to y and a dot denotes differentiation 
with respect to t. The subscript marks quantities at y = (on the brane). Also [A\ = 
A{0^) — A{0^) denotes the jump of the function A across y = 0. Assuming Z2-symmetry 
about the brane for simplicity, the junction conditions (12)-(14) can be used to compute a', 
n' and 0' on two sides of the brane. The energy-momentum conservation equation (8) on 
the brane becomes 

P& + 3— (pfc+Pb) = 00- (15) 

ao 00 

Because of the presence of the time-dependent bulk scalar field and 0-dependent couplings 
in the standard model lagrangian, the right hand side of the above equation is non-zero 
and shows the amount of energy non-conservation (due to bulk-brane energy-momentum 
transfer) of the matter fields on the brane. 

2.2 DGP braneworld cosmology with a bulk scalar field 

We use the methods presented in Refs. [16,17] (see also [9]) to obtain a special class of 
solutions for a DGP braneworld cosmology with a bulk scalar field. In this respect, following 
[17], we introduce the quantity F as a function of t and y as follows 

So, the components of the Einstein tensor can be rewritten in the following simple forms 



5 



In the presence of the bulk scalar field, the left hand sides of these two equations are not 
the same. But for special class of solutions with = 0(a), they are equivalent and in this 
case F — F{a). In this situation, both (17) and (18) then reduce to 



da) 



+ 



.a) 



+ 4pb}6{y) + 6F+-a 



3 dF 



2 da a 



^ = 0. (19) 



We choose a Gaussian normal coordinate system so that b'^{y,t) — 1. Also we assume that 
t as a proper cosmological time on the brane has scaled so that no — 1. By adopting a Z2 
symmetry across the brane, equations (13) and (16) yield the following generalization of the 
Friedmann equation for cosmological dynamics on the DGP brane 



/do\2 1 2 2fi;l k 2kI k\ 1 ^ k ^ 



(20) 



We note that here we consider only the normal, ghost-free branch of the solutions. Since we 
consider = 0(a), the field equation (2) reduces to 



^ aa^ a da'' 



{G 



Gi)- + — 



dcp dV 



da d(f) ^ y/~9 



Substituting equations Gab — K'^Tab and (19) into (21), we find 

2a\K'i„/ d(j)\2 (3Kl/d\2 3kI / k 
n\ ^a' kI 



F{a—) + 



da'' 



3 dF 

7:0'— r- 
2 da 



-h 5Cb{(p) 



0. 



(21) 



k}d(j) 
a j da 



(22) 



^ 50 

Thus the original partial differential field equations have been reduced to an ordinary differ- 
ential equation. 



2.3 Supergravity-style solutions 

In order to generate some solutions of the field equations, we introduce a special supergravity- 
style potential, V(0), as follows [9] 

Assuming k — 0, the field equations (19) and (22) are satisfied if 

F = ^W', (24) 



d0_ 3_dW 



Here, W is referred to as a superpotential, but we note that supergravity is not required for 
solutions (23)-(25) to satisfy the equations (19) and (22). Now we can rewrite the Priedmann 
equation (20) in terms of W as follows 



a2 



«0 



1 2 I ^^4 



24 



\ 



Kit: 



(26) 



From equations (25) and (26) we find the the time variation of the scalar field on the 
brane as 



1 2 , 



2«| 

4\ 



«5 



o 4Fb 3g ( 




(27) 



The jump conditions (12) and (14) are related via equation (25). Consistency of jump 
conditions for a Z2-symmetric DGP brane is guaranteed if 



^A(0) 



k\ do 



6 dW" 
4W d(f) ^ 



(28) 



So, the energy on the DGP brane is not conserved in this setup. Combining equations (15) 
and (28), one finds the energy conservation equation in terms of 1^ as follows 



pb + ^—{pb + pb) 

do 



3 dl 

Pb 2~ 

flg 



Wo 
Wo 



(29) 



We note that this equation in the case of RSII braneworld scenario has a simpler structure 
(see for instance, equation (4.7) of Ref. [9]). In our case, due to the presence of two terms 
on the right hand side of equation (29), there will be new possibilities with different cosmo- 
logical implications. To discuss the status of the conservation equation in different cases, we 
define the parameter X so that X = -% for simphcity. 



There are three possibilities as follows: 
A: pb> X 

For negative values of energy will leak off the brane. For positive values of this 
quantity, energy will flow from the bulk into the brane. 

B: pb<X 



In this case for positive values of energy will leak off the brane and for negative 



values, energy will transfer from the bulk into the brane. 
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C: pb^X 



In this case the right hand side of equation (29) vanishes and the sign of ^ is not im- 



portant, thus the energy on the brane is conserved. 



We note that for all mentioned cases, if ^ = 0, we have conservation of energy on the 
brane. Now we consider the simplest generalization of the brane energy density as 



Pb = Wop 



(30) 



where p is proportional to the energy density of the ordinary matter on the brane. This 
generalization has its origin in the fact that matter Lagrangian on the brane, that is Cb{(/)), 
depends on the bulk scalar field, 0. The effective Fricdmann equation (26) then becomes 



-i^ -kIWoP+^-^. 



4 + \4WoP+§Wl 
Kt 3 



36 



(31) 



This equation implies that if Wq tends to oo at late time, the cosmological evolution on 
the brane is not generally compatible with observations. On the other hand, for Wo — 
0, the result is not compatible with observation too. Only for a constant Wo we find 
a viable DGP-like cosmology in this simple generalization. In this case one obtains a self- 
accelerating solution which explains the late-time cosmic speed up. We note that a phantom- 
like prescription can be realized in this case in the same way as has been shown in Ref. [18] 
a canonical scalar field on the brane. 

In the next section, we will consider some specific examples of superpotential and we will 
discuss their cosmological implications. Specifically, we will consider the evolution of W in 
order to study the status of the conservation equation. 



2.4 Some specific examples 

We consider the following exponential form of the superpotential [9] 



W 



+ S- 



(32) 



ai (X2 

where s = ±1 and a\ > \a2\- For a2 — 0, we use the following form of the superpotential 



W 



+ — 

ai ^5 



(33) 



The corresponding potentials obtained from (23) are 



3al 



-2ai(f> 



+ 1 



San 



3aia2 



(34) 
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and (for 0:2 = 0) 



V 



1 - 



„ ^„ . . (35) 

Scef j 3q;i 3 J 

For \/ bounded from below, only some values of the parameters are allowed [9]. For W as 
given by (32), one can solve equation (25) to find 



ln( 



In 



— se 



-0121 



where a* is an arbitrary constant. For a2 — 0, equation (25) is solved by 



3ai 



(«:5</> + se"''^)- 



(36) 



(37) 



Wc emphasize here that although the results of this section (equations (36) and (37)) seems 
to be formally the same as the results obtained in [9] for Randall-Sundrumm II (RSII) 
braneworld, but one should remember that our Friedmann equation (26) and hence the 
Hubble parameter and scale factor differ from corresponding quantities in Ref. [9]. With 
these new quantities, equations (36) and (37) differs essentially from corresponding equations 
(5.5) and (5.6) of Ref. [9] for RSII case. 

Depending on the choice of s and the sign of 0:2, there will be a variety of cosmological 
evolution on the brane with several interesting implications. To illustrate further, in which 
follows, we will separate each of these subcases for some values of ai and a2- Then we will 
consider the evolution of scalar field versus the scale factor in normal DGP branch of the 



scenario. 



2.4.1 a2 > 0, s = +1 



In this case, as a goes from to 00, scalar field rolls down from either +00 or —00 to 0. 
In figure 1, we show the evolution of the scalar field versus the scale factor explicitly. Wc fix 
«! in a constant value and choose three different values for a2. As we see, when the value 
of a2 grows, the slope of the curve increases. This increase means that for larger a2, the 
scalar field varies faster. Also, for negative 0, variation of a2 has no significant effect on the 
slope of the curves since all curves coincide in this case. Now we look at the status of the 
continuity equation in this case. 

Jit If we consider the case that starts from —00, then as (p goes to (and a increases), 
superpotential evolves from infinity to a constant value with 1^ = (see figure 5). During 
this evolution, ^ < 0, and therefore from (29) we will have the following conditions: 

A: If ph > X, the right hand side of equation (29) becomes negative and this indicates 
that energy leaks off the brane as scale factor increases and the universe expands. This situa- 
tion continues until and W tend to zero. Then, there will be no energy leakage off the brane. 
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a 



0-2 = \ a2 = 2 a2 = 4 



Figure 1: Evolution of the scalar field with respect to the scale factor for the case 02 > 0, s = +1. 
As the scale factor evolves from to 00, scalar field rolls down from either +00 or —00 to 0. 

B: If Pi < X, the right hand side of equation (29) becomes positive and the energy is 
sucked onto the brane as universe expands. This situation continues until constant value of 
W {yV = 0) in = 0. Then, there will be no energy suction onto the brane. 

i/t On the other hand, if we consider the case that starts from +00 and tends to zero 
as a goes to infinity, the superpotential evolves from +00 to a constant value where W = 
During this stage, ^ > 0. In analogy to the previous case, there are two possibilities: 

A: If < X, the right hand side of equation (29) becomes negative. This means that 
as the universe expands, energy leaks off the brane until W = 

B: If ph > X, the right hand side of the equation of continuity becomes positive and the 
energy sucks onto the brane as scale factor increases. As soon as W vanishes, the energy 
suction stops. 

2.4.2 a2 > 0, s = -1 

As can be seen from figure 2, for this choice of parameters, evolves from +00 to —00 
or vice versa. During this evolution, the scale factor varies from zero to a maximum value 
and then reduces again to zero. For a constant value of ai, increasing the value of 02 leads 
to increasing in the slope of the curves when (p is positive. For negative 0, corresponding 
changes are not significant. Regarding to the continuity equation we have the following 
possibilities: 
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Figure 2: Evolution of tlie scalar field with respect to the scale factor for the case 02 > 0, s = — 1. 
can evolve from +00 to —00 or vice versa. During this evolution, the scale factor changes from 
zero to a maximum value and then reduces again to zero. 

^ Equation (32) shows that at first, when <^ is +00 or — cxd, W is —00 or +00 respectively 
(see also figure 5). When starts to decrease from infinity, we have ^ > 0. This situation 
continues until tends to a special value at which scale factor reaches its maximum value 
and d = 0. This implies that: 

A: For ph < X, energy leaks off the brane as scale factor increases and the universe 
expands. This leakage of energy continues until the scale factor reaches its maximum value 
and d = (see figure 2). At this point, W vanishes and there is nothing left on the brane. 

B: If ph > X, the situation is very different. When is +00 and ^ > 0, the right hand 
side of equation (29) is positive. This means that energy is sucked onto the brane whereas 
the universe expands. This suction persists until scale factor reaches its maximum value. 

^ After scale factor reaches its maximum value, the universe begins to re-collapse and 
scale factor decreases then. In this stage as and W goes to —00 and +00 respectively, 
^ < 0. This continues until the scale factor returns to zero where the superpotential tends 
to infinity. 

A: If pb < X, the right hand side of equation (29) becomes positive. This indicates that 
as the scale factor decreases and the universe re-collapses, energy is sucked onto the brane. 
This suction continues until the scale factor tends to zero where W becomes infinity. 

B: If ph > X, the energy leaks off the brane until the scale factor returns to zero where 
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the superpotential tends to infinity. 

Jl» When (f) starts to increase from — oo, we have ^ < 0. This situation continues until 
4> tends to a special value at which scale factor reaches its maximum value and d = 0. This 
implies that: 

A: For ph > X, energy leaks off the brane as scale factor increases and the universe 
expands. This leakage of energy continues until the scale factor reaches its maximum value 
and d = 0. At this point, W vanishes. 

B: li pt, < X, the right hand side of equation (29) is positive. This means that energy is 
sucked onto the brane whereas the universe expands. This suction persists until scale factor 
reaches its maximum value. 

^ After scale factor reaches its maximum value, the universe begins to re-collapse and 
scale factor decreases then. In this stage as (f) goes to +oo, ^ > 0. This continues until the 
scale factor returns to zero where the scalar field tends to infinity. So: 

A: If pb > X, the right hand side of equation (29) becomes positive. This implies that 
as the scale factor decreases and the universe re-collapses, energy is sucked onto the brane. 
This suction continues until scale factor tends to zero where (p becomes infinity. 

B: If ph < X, the energy leaks off the brane until the scale factor returns to zero where 
the scalar field tends to infinity. 

2.4.3 q;2 < 0, s = -1 

In this case, the scalar field evolves from — oo to +oo, as the scale factor goes from zero to 
oo. Figure 3 shows this behavior. One can see from this figure that by decreasing the value 
of a2, the slope of curves increases for positive scalar field. However, for negative scalar field 
this change is not significant and there is no considerable shift in curves. Regarding to the 
continuity equation, the following issues are in order: 

^ At first, when the scalar field is — oo, the superpotential is +oo. Then, by increasing 
the values of scale factor and scalar field, the superpotential decreases to zero as (p ^ oo. 
During this evolution, we have ^ < 0. 

A: If pb > X, as scale factor increases and the univers expands, energy leaks off the 
brane and this leakage continues until W tends to zero at late time. 

B: li pb < X , from equation (29) we find that as the universe expands, the energy is 
sucked onto the brane. 
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Figure 3: Evolution of tlie scalar field with respect to the scale factor for the case with 02 < 
and s = — 1. As scale factor goes from zero to 00, the scalar field evolves from —00 to +00. 

2.4.4 as < 0, s = +1 

For this choice of parameters, the solution starts from = —00 and a = 0. The scalar 
field increases with scale factor until some value of a where d = 0. At this point a reaches a 
relative maximum. After that, the scale factor begins to decrease and scalar field continues 
its growing until it tends to zero. At this point scale factor reaches its relative minimum, 
where a is zero. Then, grows again with scale factor towards infinity. This behavior can 
be seen in figure 4. Further, this figure shows also that decreasing the value of 0^2 leads to 
increasing in the slope of curves for positive 0. For negative 0, change in the values that 
a2 attains has no significant effect on the evolution of 0. Once again, regarding to energy 
conservation in this case we arrive at the following points: 

J(t Firstly, when is —00, the superpotential is +00. As increases and a goes to its 
relative maximum, W tends to zero. During this stage, ^ < 0. So we have: 

A: If ph > X, the right hand side of conservation equation (29) becomes negative and 
energy leaks off the brane with expansion of the universe. This feature continues until W 
vanishes at the relative maximum of the scale factor. 

B: If pb < X, the right hand side of equation (29) becomes positive and energy is sucked 
onto the brane as scale factor increases. This situation continues until the scale factor reaches 
its relative maximum. 

J(> Secondly, when a decreases to its relative minimum and tends to zero, the superpo- 
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Figure 4: Evolution of the scalar field with respect to the scale factor for the case with 02 > 
and s = +1. The scalar field increases with scale factor until a reaches a relative maximum. After 
that, the scale factor begins to decrease and the scalar field continues its growing until it tends to 
zero. Then, grows to infinity when a tends to infinity. 

tential reaches its minimum (here it is a negative value) and W vanishes. During this stage, 
f > 0. We find that: 

A: If pf, > X, the right hand side of equation (29) becomes positive. As the scale factor 
decreases and the universe re-collapses, energy is sucked onto the brane until = 0. 

If < X, the right hand side of conservation equation becomes negative. So, the 
energy leaks off the brane and this leakage persists until the minimum value of superpotential 
is achieved. 

i/t Thirdly, when a starts to increase from its relative minimum towards infinity, the su- 
perpotential changes from its minimum to zero when goes to +00 . During this evolution, 
we have ^ < 0. Here also we have: 

A: If pfe > X, the right hand side of equation (29) becomes negative and as a and tend 
to infinity, the energy leaks off the brane. 

B: If pfe < X, the right hand side of equation (29) becomes positive. Here, as universe 
expands, the energy is sucked onto the brane. 
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Figure 5: Evolution of the superpotential versus the scalar field. This evolution plays important 
role in the status of the energy conservation. 

2.4.5 = 0, s = ±1 

In these cases, for s = +1, scalar field evolves from — oo to +oo, while for s = — 1 it 
evolves from — oo to +oo or vice versa. The value of s determines status of this evolution. 

Jlk If s = +1, (f) starts from — oo at a = and continues to +oo as a goes to infinity. When 
(j) is — oo, the superpotential is +oo. As scalar field increases, the superpotential decreases 
to a constant value when cj) tends to infinity. During this stage ^ < 0. So we have: 

A: If pb > X, as a and go to infinity, the energy leaks off the brane. This leakage stops 
when W reaches a constant value. 

B: If ph < X, the right hand side of equation (29) becomes positive. Here, as universe 
expands, the energy is sucked onto the brane. 

J|k If s = — 1, starts from — oo at a = and increases until a reaches its maximum. At 
this point d and W tend to zero. In this situation ^ < and we have: 

A: If ph > X, as a and (f) go to infinity the energy leaks off the brane. This leakage 
continues until W tends to zero. 

B: If pf, < X, as universe expands, the energy is sucked onto the brane. 

Just after d becomes zero, the scale factor starts to decrease and the universe re-collapses. 
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Figure 6: Evolution of scalar field respect to scale factor. If s = —1, scale factor has a maximum. 
But for s = +1 it continues to infinity. 

When a goes from maximum to zero, the scalar field continues its growth and tends to +oo. 
During this stage, W decreases and we have ^ > 0. So the following subcases are in order: 

A: If ph < X, as the universe re-collapses, the energy is leaks off the brane until W tends 
to a constant value. 

B: li ph > X a.s a goes to zero, the energy is sucked onto the brane. This suction persists 
until the scalar field reaches +oo. 

^ For s = — 1, there is another situation that starts from +oo at a = and decrease 
until a reaches its maximum. At this point a and W tend to zero. In this situation ^ > 
and we have: 

A: li Ph < X, as a increase and decrease, the energy leaks off the brane. This leakage 
continues until a reaches its maximum value. 

B: If pfc > X, as universe expands, the energy is sucked onto the brane. 

Just after a becomes zero, the scale factor starts to decrease and the universe re-collapses. 
When a goes from maximum to zero, the scalar field continues its decrease and tends to —oo. 
During this stage, W increases and we have ^ < 0. So the following subcases are in order: 

A: If pb > X, as the universe re-collapses, the energy is leaks off the brane until W tends 
to a constant value. 
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B: If pf, < X as a goes to zero, the energy is sucked onto the brane. This suction persists 
until the scalar field reaches +00. 




Figure 7: Evolution of the superpotential versus the scalar field for 02 = 0. If s = —1, the 
superpotential will pass through a zero and then tends to a constant . If s = +1, the superpotential 
has no root. 

In summary, assuming a simplest generalization of the brane energy density {ph = Wop), 
implies that if Wq tends to +00 at late time, the cosmological evolution on the brane is 
not generally compatible with observations. On the other hand, for Wq = 0, the result is 
not compatible with observation too. Only for a constant Wq we find a viable DGP-like 
cosmology in this simple generalization. So, we found in this analysis that the case with 
^2 > and s = +1 gives a viable DGP-like cosmology. We should note also that in this 
case, the situation that the energy leaks of the brane during universe expansion is more 
reliable. 

After a detailed discussion on the energy conservation and possible cosmological dynamics 
on the brane, in the next subsection we determine cosmological dynamics in the bulk. 



2.5 Bulk solutions 

By solving the non-zero off-diagonal components of the Einstein field equations, we are able 
to determine the y-dependence of the metric or the scalar field (p in the bulk. To do this 
end, we assume the fifth dimension is static. In other words, 6 = and so we can adopt the 
gauge 6 = 1. When = 0(a), we have 

(n'a ha' a'\ ^ ^ , . ( d(t)V 

\ na ba a \ da 
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Then n can be expressed in terms of a according to the following relation 



^ = /3(t)^exp 
n 6 



/ a( ^ 1 da 

da 



(39) 



where ^ is a function of time alone and has no dependence on y. liV has a supergravity-like 
form as (23), we can rewrite the above expression as 



(40) 



We set no = 1, then (26) imphes that 



/3 = ^ 



724 UpkKl 



72Ki 



A 



4 + Up- + 5w? 

Kt 3 



nl/2 



36 



(41) 



To have a simple comparison, we note that in the RSII case this quantity is given by f]^^^^^"* 



0'0\l pI/Wq — 1 (see [9] for instance ). Inserting equation (40) into (16) leads to the following 
differential equation 

(42) 



do 



We take ai = q;2 = so that the superpotential used in subsection 2.4 can be simplified 
to 

\ 2 

a* 



W — — cosh (q;i0) — 



+ 4 



a 



(43) 



Therefore the differential equation has a general solution of the form ( see [17] for instance ) 

a^ — A cosh jiy + B sinh jiy + C, (44) 

where 



2CK5 



(45) 



and the coefficients A, B and C are functions of time. 

The Z2 symmetry across the brane imposes the relations Aj^ = A_ = A and B^ = B_ = 
B between these coefficients on the two sides of the brane [17]. Prom equations (12) and 
(13) we can find the following relations between coefficients 



B/i 
A + C 



3 AC^ Ojq 



(46) 



and 



2^ 



Bp 



(5kIW^^2. Dp 



A + C 



{A + cy 



^(3pfe + 2p6) + ^ 



a. 



_ zo _2: 
al ao 



aQfiQ 



+ 2 



ao 



(47) 
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Using equations (26), (40) and (44) and taking the limit y — 0, one can determine the 
coefficients as follows 



C 




ao 



18 KiXaoJ 3 



(48) 



2a| 



B 



+ 



do 
ao 



(49) 



(50) 



Substituting these coefficients into the general solution (44), we obtain the following expres- 
sion for the bulk behavior of the scale factor 



a — an 




4 2 / • \ 2 2 



As we mentioned previously, W ca and therefore, the singularities in this model occur 
if — >■ oo. In figure 8, we plotted the scale factor versus and y. This figure shows that 
for some values of pf, and y, there are points in the parameter space that the scale factor 
vanishes. The position of these singularities varies with cosmic time on the brane. At these 
points, 0', (f) and W all are singular and is so Tm- Thus the singularities at = are 
naked curvature singularities. Nevertheless, if we have a compact bulk or more other branes 
at suitable distances from our brane, we can avoid these singularities. Also, if we have the 
solutions with bounded W ^ it is possible to elusion from the curvature singularities. Note 
that in ploting figure 8, we have set Wq = 2ph. 

If we take ai = —a^ = the differential equation (42) has a general solution of the 



form 



Ay^ + By + C. 



(52) 



Following the same procedure as above, one can find the following expression for the bulk 
scale factor 



1 - 




2 k"^ (<-2 





(53) 



Mo J 18 ° nlxaQ, 

Figure 9 shows variation of a^ versus and y. As this figure shows, the bulk scale 
factor vanishes in some points. Since Tab is also divergent in these points, there are naked 
curvature singularities in the bulk. 
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Figure 8: Evolution of the bulk scale factor. 




Figure 9: Evolution of the scalar field. For some values of pi, and y, it becomes zero. In this points 
there are naked singularities. 
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3 Bulk phantom scalar field 



3.1 Field equations 



In this section we consider a DGP-inspired braneworld model that the bulk contains a phan- 
tom scalar field. For a phantom scalar field, the sign of the kinetic energy term is opposite 
of the canonical scalar field case studied in the previous section. The five-dimensional action 
for this model can be written as follows 



(54) 

Variation of the above action with respect to the phantom scalar field and also the metric, 
gives the following equations of motion. These equations differ from equations (2) and (3) 
for canonical scalar field in the sign of the scalar field dependent terms 



d(j) 



'—h dCb{(j)) 
^ d(t> 



(55) 



2rpA{hrane) 



(56) 



In the presence of the phantom scalar field, the Gauss-Codacci junction conditions are 
the same as for the canonical scalar field case. Also, if we assume the metric to be as (9), the 
jump conditions and the energy-momentum conservation equation, are the same as given in 
section 2.1. 

In order to find some solutions of the five dimensional field equations in this case, we 
use the quantity F and the component of the Einstein tensor used in section 2.2. By using 
equation (16), for special class of solutions with (p = 0(fl) and so F = F(a), equations (17) 
and (18) reduce to 

2^,. ,s i^T^, d(p,2 f3K?/d\2 3/t? / A; \ r. ] ^, . „^ 3 dF 3k ^ . 



Also the field equation (55) reduces to 



^ aa"^ a da^ 



+ 



da 



^ ^A(0) 
5(1) 



(58) 



By using Gab — k^Tab and equation (58), equation (57) gives 



F(a 



d , 
d^' 



2a 



W(-) +^ 

K/^ Qj K/^ 



-) + 4pb\6{y) + 6F+-a 



da 



= 0. 



1 d4> 
a j da 



(59) 



This is the effective field equation for a bulk phantom scalar field in the DGP setup. In 
which follows, we present some supergravity-style solutions for this effective field equation. 
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3.2 Supergravity-style solutions 

In this subsection, to present some solutions of the field equations, we use a supergravity- 
style potential introduced earlier in section 3.2. Assuming /c = 0, the field equations (57) 
and (59) are satisfied if 

(60) 
(61) 

(62) 



36 ' 



3 dW 



a— — 



da K^W d(j) ' 

respectively. The Priedmann equation on the brane now takes the following form 



"O _ ^ 2 



2k| 



A 



4 1 2 



Ave; 



Also, time variation of the phantom scalar field now is given as follows 
'1 



24 

4 



24 

Kit: 



A 



9 



dWp 
~d^ 



(63) 



We note that the energy conservation equation is the same as given by equation (29). 



3.3 A phantom field superpotential 

As discussed in Ref. [19], the curvature of the universe grows toward infinity within a finite 
time in the universe dominated by a phantom fiuid. In the case of a phantom scalar field, 
this Big Rip singularity may be avoided if the potential has a maximum. Here we consider 
the following form of the superpotential which has the mentioned property 



W{(t)) = WA cosh 



, acf) , 



(64) 



This superpotential has a maximum at = and tends to zero when the phantom scalar 
field grows to infinity (see figure 8). The corresponding potential obtained from equation 
(23) is 



V{<P) 



1 / a 



4 ^nipi 



(^tanh(— )J - 



rripi 



cosh 




The solution of equation (61) when W is given as (64) is as follows 



a(f) 
^rupi 



(65) 



(66) 



where a* is an arbitrary constant. This relation shows that for a bulk phantom scalar field, 
the result of the superpotential approach for evolution of the scalar field versus the scale 
factor, is different for even and odd values of a. For odd values of a, scalar field varies from 
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Figure 10: Evolution of the phantomic superpotential (64) with respect to the phantom scalar 
field. This superpotential has a maximum at (f) = 0. 




K. = 1 , m , = 1 
5 pi 



■ a =2 a =4 ~ ■ ~ a =6 



Figure 11: Evolution of the phantom scalar field versus the scale factor for even values of a. 
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Figure 12: Evolution of the phantom scalar field with respect to the scale factor for odd values of 
a. 

zero to -f oo as scale factor grows to infinity. But, for even values of a, as scale factor goes 
to infinity, scalar field varies from zero to both ±oo. 

As figure 8 shows, W has its maximum value at = 0. In this point, since Wq = too, 
the right hand side of equation (29) vanishes and we have the energy conservation on the 
brane. 

Jft When goes to +oo, W decreases until it tends to zero when scalar field reaches 
infinity. During this course, ^ < 0. This implies that: 

A: For > X, energy leaks off the brane as scale factor increases and the universe ex- 
pands. This leakage of energy continues until the scale factor reaches infinity and W tends 
to zero. After that there is no leakage of energy-momentum off the brane. 

B: However, if < X, the situation is different. Since ^ < 0, so the right hand side 
of the conservation equation is positive. This means that energy is sucked onto the brane 
whereas the universe expands. This suction persists until scale factor and the scalar field 
reach infinity. After that the suction of energy onto the brane will be stopped. 

J(t As (j) goes from to — oo, W decreases until it tends to zero when scalar field reaches 
infinity. During this course, ^ > 0. This implies that: 

A: For pf, < X, energy leaks off the brane as scale factor increases and the universe 
expands. This leakage of energy continues until W tends to zero. After that there is no 
leakage of energy-momentum off the brane. 
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B: However, if pi, > X , the right hand side of the conservation equation is positive. This 
means that energy is sucked onto the brane whereas the universe expands. This suction 
persists until scale factor and the scalar field reach infinity. After that the suction of energy 
onto the brane will be stopped. 

We note that since in this case W tends to zero at late time, the model with bulk scalar 
field with superpotential as given by equation (64), is not a viable DGP-like cosmology. 



3.4 Bulk solutions with phantom field 



To find the y dependence of the metric (or (6), we proceed the method used in section 5.2. 
As before, we assume the fifth dimension is static (6 = 0). When = 4>{a), we have 



n a 



ha' 



a' 



05 



na'^ ha a 



5-^05 



2 / • I 

—K^a a 



da 



So, we can find the expression of n in terms of a 



a 
n 



= ^(t)^cxp 
6 



da 



da 



(67) 



(68) 



As before, /3 has no y-dependence. If V has a supergravity-like form as (23), one can simplify 
the above equation to 



n 6 



For riQ — l, equation (62) gives 

"12^1^6 



Wo 



+ 



By substituting (69) into equation (16), we find 



3 36 ° 



(69) 



(70) 



(71) 



Following the same procedure as adopted in subsection 2.5, and by redefinition of W as 
equation (64), we find the following expression for the bulk behavior of the scale factor 



+«0 



1 - 



p 



1^5 , 1^5 f^O 
3 ^4 V Oo 



sinh {p\y\ 



(72) 



We note that for a bulk phantom scalar field with a bounded superpotential, as given by 
(64), the model has no singularity. 
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4 Summary and Discussion 



In this paper, we have studied cosmological dynamics in a DGP setup with a bulk canoni- 
cal/phantom scalar field. Bulk scalar field is motivated in several context: to stabifize the 
distance between two branes in the Randall-Sundrum two-brane model, a way to solve the 
cosmological constant problem, and inflation driven by bulk scalar field without inflaton on 
the branc. Generally, the scalar field living in the bulk affects the cosmological dynamics on 
the brane considerably. The evolution of this field has interesting cosmological effects and 
can give rise self-acceleration and phantom-like phase even in the normal DGP branch of 
the model in some appropriate situations. In this paper, we have generalized the work of 
Ref. [9] to the DGP setup. We considered an extension of the DGP scenario that the bulk 
is non-empty and contains a canonical or phantom scalar field. Since the self-accelerating 
DGP branch has ghost instabilities, we restricted our study to the normal DGP branch of 
the model. The bulk equations of motion are derived and some special classes of the solu- 
tions are presented. We determined also the evolution of the brane when the potential of the 
scalar field takes a supergravity-like form. Some clarifying examples along with numerical 
analysis of the model parameters space are presented in each step. The importance of this 
work is that bulk scalar field in the DGP setup was not studied in supergravity-style analy- 
sis. Also our detailed study in this framework has revealed some yet unexplored aspects of 
cosmological dynamics of the bulk scalar field in DGP setup. We have extended this study 
to the case that the bulk contains a phantom scalar field too. Our strategy to perform the 
mentioned analysis was as follows: 

First of all, from a five-dimensional action for a DGP-inspired braneworld model with 
a bulk canonical scalar field, we found the bulk equations of motion, jump conditions and 
the Gauss-Codacci equations. Then, using these equations, we achieved the effective energy- 
momentum conservation equation on the brane. We saw that because of the presence of 
the time-dependent bulk scalar field and 0-dependent couplings in the standard model La- 
grangian, the right hand side of the continuity equation is non-zero, showing the amount of 
energy non-conservation (due to bulk-brane energy-momentum transfer) of the matter fields 
on the brane. To obtain a special class of solutions for a DGP braneworld cosmology with 
a bulk scalar field, we used the methods presented in Refs. [9,16,17]. Using that method 
and introducing the quantity F as a function of t and we reduced the original partial 
differential field equations to an ordinary differential equation. In order to generate some 
solutions of the field equations, we introduced a special supergravity-style potential V"(0), in- 
cluding the superpotential W . With this supergravity-style potential, we derived the energy 
conservation equation on the brane in terms of W . In our model, due to the presence of two 
terms on the right hand side of the conservation equation (29), there was new possibilities 
with different cosmological implications. If the right hand side of the conservation equation 
becomes negative, the energy leaks off the brane. However, if the right hand side of this 
equation becomes positive, there is energy suction onto the brane. For vanishing right hand 
side of the conservation equation, energy is conserved on the brane. We considered some 
specific examples of superpotential and discussed their cosmological implications. Evolution 
of the scalar filed versus the scale factor, and evolution of W in terms of the scalar field are 
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discussed fully to study the status of the conservation equation on the brane. Assuming the 
simplest generalization of the brane energy density {ph — Wqp), implies that if Wq tends to 
+00 at late time, the cosmological evolution on the brane is not generally compatible with 
observations. On the other hand, for Wq = 0, the result is not compatible with observation 
too. Only for a constant Wq we find a viable DGP-hke cosmology in this simple general- 
ization. So, we found in those examples that the case with a2 > and s = +1 (subsection 
2.4.1) is a viable DGP-like cosmology. Of course in the mentioned case, the situation that the 
energy leaks of the brane during universe expansion is more favorable. We also determined 
the bulk behavior of the metric. We saw that since W is not bounded, the naked curvature 
singularities in the bulk are present. However, with a compact bulk or existence of more 
other branes at suitable distances from our brane, we can avoid these singularities. Also, 
for the solutions with bounded W, it is possible to elusion from the curvature singularities. 
We continued our treatment by considering the bulk phantom scalar field. From the five- 
dimensional action for a DGP-inspired braneworld model with a bulk phantom scalar field, 
we found the bulk equations of motion. In the presence of the bulk phantom scalar field, 
the Gauss-Codacci junction conditions and also the energy-momentum conservation equa- 
tion, are the same as for the canonical scalar field case. Introducing quantity F in terms 
of metric components as before, we achieved the effective field equation for a bulk phantom 
scalar field in the DGP setup. To present some solutions of the field equations, we used a 
supergravity-style potential introduced earlier in the case of canonical scalar field but with 
a new superpotential. We considered the evolution of the scalar field versus the scale factor 
and the evolution of W versus (j) in order to study the status of the conservation equation on 
the brane. Since in this case W tends to zero at late time, this case is not a viable DGP-like 
cosmology. Also, we determined the bulk behavior of the metric. Since in this case W is 
bounded, there is no singularity in the brane. 

Finally the following issues are important to note: 

• To be a cosmologically viable scenario, this model with bulk scalar field should explain 
at least the late-time cosmic speed-up on the brane. As we have shown, this model 
accounts for cosmic speed-up on the brane in some specific situations. For instance, 
the case with a2 > and s = +1 in subsection 2.4.1 gives a viable DGP-like cosmol- 
ogy explaining late-time cosmic acceleration. Also it is possible to realize an effective 
phantom-like prescription on the brane without need to phantom matter in the same 
way as has been done in Ref. [18] for a canonical scalar field on the brane. 

• As has been mentioned by Flanagan et al. (in Ref. [8]), the principal problem as- 
sociated with the introduction of a bulk scalar field in models such as the present 
one is that generic stationary solutions contain timelike curvature singularities in the 
bulk at finite distances from the branes. One approach to overcome this problem is 
to simply orbifold or otherwise compactify the fifth dimension in so that the singu- 
larity is never encountered. A second approach is to carefully choose the scalar field 
potential in such a way that the occurrence of singularities is prevented. Here we have 
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adopted the second strategy by choosing bounded superpotentials. For instance, with a 
bulk phantom scalar field with a bounded superpotential, the model has no singularity. 

• As has been shown in Refs. [12] and [13], in the presence of a bulk scalar field, real- 
ization of the inflation is possible without inflaton field on the brane. In fact, inflation 
can be driven just by a bulk scalar field. In our setup, it is possible to realize infiation 
in the same way as has been adopted in Refs. [12] and [13]. In fact, the late-time 
behavior of the bulk scalar field can be treated by analyzing the property of a retarded 
Green function. Including the lowest order back-reaction to the geometry, this late- 
time behavior can be well approximated by an effective four dimensional scalar field on 
the brane. As has been shown by Himemoto et al. in Refs. [12,13], the mapping to the 
four- dimensional effective theory is given by a simple scaling of the potential with a 
re-definition of the field. This effective four-dimensional field can drive infiation on the 
brane. This issue is under study and will be addressed in one of our forthcoming report. 
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